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spacetime. Due to the super-radiant effect of the background geometry, the bound 
state of the scalar field is unstable, and its amplitude grows in time. By solving the 
Klein-Gordon equation of the scalar field as an eigenvalue problem, we numerically 
obtain the growth rate of the amplitude of the scalar field. Although the dependence 
of the scalar field mass and the scalar field charge on this growth rate agrees with the 
result of the analytic approximation, the maximum value of the growth rate is three 
times larger than that of the analytic approximation. We also discuss the effect of the 
electric charge on the instability of the scalar field. 
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§1. INTRODUCTION 



The propagation of waves in black hole spacetime has been studied by many researchers 
in the context of black hole perturbations P Super radiance is one type of interesting 
phenomena related to this topic PEM3J3 Considering the situation that a wave impinges on 
a black hole, the incident wave is partially reflected by the potential barrier of the black hole, 
due to the centrifugal force, and is scattered to infinity, while part of the wave penetrates 
the potential barrier and is absorbed by the black hole. Thus, one would expect that the 
amplitude of the scattered waves is always smaller than that of the incident wave. However, 
this is not necessarily true in the Kerr-Newman geometry. If the frequency of the incident 
wave satisfies the so-called super-radiant condition, the reflected wave is amplified, and its 
amplitude becomes larger than that of the incident wave; that is, the reflection rate can be 
greater than unity. Using this amplification mechanism, it is possible to extract the rotation 
and the electro-static energy of Kerr-Newman black holes. 

As an application of the super radiance, Press and Teukolsky proposed the "black hole 
bomb" P They considered the situation in which mirrors surround a Kerr black hole. Then, 
the scattered wave with super-radiant amplification is reflected back to the black hole by 
the mirrors, and the wave is amplified again. In this way, the amplitude of the wave grows 
exponentially in time and becomes unstable. Damour et a/P have shown that a black hole 
bomb can be realized by using a charged massive scalar field in Kerr-Newman spacetime. 
The effective potential of the massive scalar field has a local minimum, and the scalar field 
wave can be confined in this potential well. Hence, the mass of the scalar field plays the 
role of the mirrors, which is necessary to cause the instability of the scalar field. According 
to their analysis, the instability of the scalar field is realized if the bound state and the 
super-radiant condition are both satisfied. 

To evaluate the growth rate of the scalar field, we must solve the Klein-Gordon equation 
in the black hole geometry. By imposing ingoing boundary conditions at the black hole 
horizon and regular boundary conditions at infinity, the problem of obtaining the unstable 
mode reduces to an eigenvalue problem. For Kerr spacetime, the growth rate of the unstable 
modes for a scalar field with a small mass /i < 1/M (where M is the mass of the black 
hole) was obtained by using the asymptotic matching method P and that for a large mass 
/i ^> 1/M was obtained by using the WKB approximation P However, these analyses do 
not cover the parameter region in which \xM ~ 1, where the growth rate is expected to take 
its maximum value. 

In this paper, we aim to obtain the growth rate of a charged massive scalar field with 
mass satisfying fiM < 1 in Kerr-Newman spacetime by solving the eigenvalue problem 
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numerically. The paper is organized as follows. In §21 we estimate the growth rate of a 
scalar field with fiM <C 1 and qQ <C 1 in Kerr-Newman spacetime using the method of 
DetweilerP Then, in JJJ we introduce our numerical method and obtain the growth rate 
caused by the instability and its parameter dependence for fiM < 1. Section 0] is devoted to 
a summary and conclusion. 

We use the units in which G = c = h = 1 throughout the paper. 

§2. Analytic approach 

In this section, we calculate the growth rate of the scalar field using the asymptotic 
matching method used by Detweiler^ in Kerr-Newman spacetime. 
The Kerr-Newman metric in the Boyer-Lindquist coordinates is 

dg2 = _ / _ ™r-C?\ df _ (2Mr-Q^e dm 

+ |*J + + (> + J + (Wr-Q>W«j ^ ^ (J) 
A = r 2 -2Mr + a 2 + Q 2 , S = r 2 + a 2 cos 2 9, 

where M is the mass, a is the angular momentum and Q is the electric charge of the black 
hole. The locations of the horizons r± are given by the roots of the equation A = 0. The 
Klein-Gordon equation for a charged scalar field with mass \x is 

( V a - iqA a ) (V a - iqA a )V = fi 2 V, (2) 
A a = I -— ,0,0, — sm 9 

where V a is the covariant derivative in the Kerr-Newman geometry, and q is the charge of 
the scalar field. Equation (j2J) is separable in terms of the spheroidal harmonics S(9): 

$ = ip(r)S(e)exp(i(-ut + m<f>)). (3) 
The radial function ip(r) satisfies the relation 

-Z\(/i 2 r 2 + A) + { (r 2 + a 2 )uj -ma- qQr} 2 ip = 0, (4) 
with 

A = /(/ + 1) - 2mau + (au) 2 + 0(a 2 fi 2 (l - uj 2 / fx 2 )), 



dr dr 
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where l,m are integers and \m\ < I. We assume I > 1. (For I = 0, there is no centrifugal 
force, and the bound state of the scalar field does not exist.) We are interested in the eigen- 
mode whose frequency is nearly equal to the mass of the scalar field, i.e. for which u ~ \i. 
As in this case we have |1 — uj 2 / fi 2 \ <C 1, we can realize the relation 0(a 2 fi 2 (l — uj 2 / /i 2 )) ~ 0, 
and then the separation constant is given by 

A = /(/ + 1) - 2mau + (auj) 2 . (5) 

We solve Eq.(@]) using regular boundary conditions at infinity and ingoing boundary con- 
ditions at the black hole horizon. To apply the asymptotic matching method, we need to 
assume that the parameters satisfy the conditions 

0(\uM\) = O(fiM) = O(qQ) = 0(e), e < 1, (6) 

and that they can thus be treated as small parameters. The angular momentum of the black 
hole is assumed to satisfy 0(a/M) = 1. 

For large r, i.e. r ^> r + , Eq. (jlj) reduces to 



d 



2 



dr 2 



2 2 2(2Mlu 2 - Mfi 2 - qQuj) 1(1 + 1) + e 
to — fi H — 



(A 1 ' 2 ^) = 0. (7) 



To satisfy the regular boundary conditions at infinity, the phase of u is required to satisfy 

< arg \/oj 2 - n 2 < it. (8) 

The solution of this equation that is regular at infinity is 

rifj = {-2ikr) l+1 e ikr U{l + 1 - v + e 2 , 21 + 2 + 2e 2 , -2ikr), (9) 

2Mu 2 - Mix 2 - qQuj 



—ik 



where U is one of the confluent hypergeometric functions The asymptotic behavior of 
the solution (jHJ) for \kr\ 1 takes the form 

(-2ikr) l+1 n . , 

rib i (10) 

V sin [tt(2Z + 2 + 2e 2 )] K ' 

X { r(-l -v- e 2 )r(2l + 2 + 2e 2 ) + " ' 

\-2ikr)- 21 - 1 _ v+l + e2)r( _ 2/ _ 2e2) + ■ ■ 
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In the region near the black hole horizon (r <C l/fi) , Eq. (j3J) reduces to 



d d 

z l~ z l~^ + 
dz dz 



P 2 - 1(1 



V> = o, (n) 



r — r + (ri + a 2 )u; — ma — qQr + 
z = , P = , 

and the solution with ingoing boundary conditions at the horizon is given by 

ip = z iP (l - z) l+1 F(l + 1, Z + 1 + 2iP, 1 + 2iP, z), (12) 

where F is the Gauss hypergeometric function. The asymptotic form of the solution (J 121) 
for 1 — z <C 1 is given by 

r(i + *p)r(2i + i) / r V 

^ r(Z + l)r(l + 1 + 2zP) \r + -r_J 

r(i + 2iP)r(-2i-i) f r y 1 - 1 

r(-l)r(-l + 2iP) \r + -r_J 1 J 

We get the asymptotic behavior described by (fTUj) and ()13|) . These expressions are valid in 
the ranges r + Cr <C 1/1; and r + < r < respectively. Consequently, for w ~ /i, we 
have the overlap region in which both the outer expansion (fTUJ) and the inner expansion (|13J1 
hold: 

r + < r < ■ (14) 

2y A* ~ ^ 

In this region, we can match the leading-order terms of the solutions (jlOj) and ()13|). and this 
matching yields 

rH-,-^r(2/ + 2) iP-htW , lpa /(* + i)r(-2i-i) 

(15) 

For uj ~ « and Zc(r + — r_) <C 1, the right-hand side of the this equation is 0((k(r + — r^)) 2l+l ). 
Therefore, we define as the value of v for which the right-hand side equals zero and write 
v = + Sv. 

For i/ = the matching condition (JTSJl yields 

r{-l - ^ - e 2 )r(2l + 2) 

r(z - + 1 + e 2 )r(-2i - 2e 2 ) 1 J 

Thus, using the property of the gamma function l/r(—n) = 0, we have 

Z - + 1 + e 2 = -n, (17) 
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where n = 0, 1, 2, • • • . From the definition of v and Eq. we find 

2Mcj(°) 2 - M/i 2 - gQcu (0) 



/(°) 



/Cj(°) 2 - (U 2 



Mfi-qQ 

I + 71 + I 



~ l + n+1, 



1/2 



1 (M\i-qQ 



2 \ l + n + 1 



fc<°> 



'a; 



(0)2 



Mn-qQ 

/ + n + 1 



(19) 
(20) 



From Eqs. (jSJ) and (|2*D*jl. to satisfy the regular boundary conditions at infinity, we must require 
the condition 

Mil > qQ. (21) 

With this condition, the effective potential has a well (see Fig.|2J), and this condition guar- 
antees that the wave is confined in this well. 

Next, we obtain 5v perturbatively. The left-hand side of Eq. becomes 



r(-l - i/°> - e 2 )r(2l + 2) 
r(/-z/(°) + l + e 2 )r(-2/-2e 2 ) 



[1 + {${1 - z/ (0) + e 2 ) - ^(-Z - z/ (0) - e 2 )}5z/] , (22) 



with 



Mz) = r\z)/r{z). 

Then, using the property of the gamma function 

Hm M = (-irv r( - n) 



4) n ~ m m! 



we find 5u as 



5i/ = 2iP {0) 



where 



2{M\i - qQ) 
l + n + 1 



-i 2/+1 



/i(r + - r_) 



x 



r(—m) 



{2l + n + 1)! 



n! 



2! 



(2/)!(2/ + r 



(23) 



n2 ri(j 2 +4P^ 2 ), 

(24) 



p(0) 



(ri + a 2 ) a/ ) — ma — qQr^ 



r + — r_ 

Here, the relation between 5v and 5u is given by 

Mil - qQ u^5u (l + n + l) 3 5u 



8v = \i- 



% A;(°) 3 " (Mp-qQ) 2 p, 



(25) 
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Thus, the real part cr and the imaginary part 7 of the leading-order contribution to the 
eigenvalue u are 



a — fi 



1 - - 



1 / Mfi — qQ x 



2 V 1+ 1 + n 



H[l- 0(e 2 )] , (26) 



1 (/ + 1 + n) 6 

If the imaginary part 7 of u is positive, the mode in question is unstable, and 7 represents 
the growth rate of the scalar field. From Eqs. (|24|) and (J27j) . the condition for the existence 
of an instability is 

p(°) = _ ( m f2 H + q$ H - a) > 0, (28) 

where the surface gravity k, the angular velocity f2 H , and the electric potential <E> H for the 
Kerr-Newman black hole are introduced as follows: 



K= l( r -±Zl=.), = *» = -gV (29) 

2 \ a 2 + r± J a 2 + ri a 2 + 



The condition of instability (|28|) coincides with that of super radiance. If the super-radiant 
condition is compatible with the condition for the existence of bound state (J21j) . the scalar 
field becomes unstable. These features of the unstable mode are consistent with the results 
of the analysis carried out by Damour et al P The most unstable mode corresponds to 
I = m = 1 and n = 0, and the value of 7 for this mode is given by 

7 = ^|M/i - gQ| 5 (a 2 + r\f{n H + q<P H - ^ + {Q H + q<P H - /i) 2 ). (30) 
For the extreme case in which a 2 + Q 2 = M 2 , 7 becomes 

7 = ^(a 2 + M 2 ) 3 |M/i - qQ\\Q H + q<P H - /i) 3 . (31) 
For a Kerr black hole (Q — 0), 7 reduces to 

7 M = — — — ^, (32) 
M 24 V ; 

which is the formula derived by Detweiler P 

The mass and charge dependences of 7 are displayed in Fig.^ The maximum value is 
7 s; 3 x 10~ 8 , realized at /iM w 0.3, ~ —0.08. The value of 7 is positive in the region 
where both the super-radiant condition > and the bound state condition M/j, > qQ are 
satisfied. Let us consider the region satisfying fiM <C 1, \qQ\ <C 1, where our approximation 
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Fig. 1. The dependence of the mass /z and the charge q of the scalar field on the growth rate 7. 
The parameter values here are a = 0.98M, Q = 0.01A/, I = m = 1, and n = 0. The solid line 
corresponds to = 0, and the dotted line corresponds to M/j, — qQ = 0. 



is expected to be good. In this region, retaining only the leading-order terms of the small 
parameters, 7 becomes 



7 



24 



\Mfi-qQ\ 5 a(M 2 - Q 2 



(33) 



and the (/x, q) dependence of 7 is determined by the factor /i 4 |M/i — qQ\ 5 . For a fixed value 
of q, 7 is an increasing function of fi, for a fixed value of fi, 7 is a decreasing function of qQ. 
Equation (J33|) shows that the main effect of the scalar field charge q is to change the depth of 
the potential well that is necessary to bound the scalar field. As the charge q increases, the 
depth of the well of the scalar field effective potential decreases, and the scalar field becomes 
less bounded. Thus, an increase of the charge decreases the growth rate of the scalar field. 

For a positive charge, i.e. q > 0, can be positive for a negative azimuthal quantum 
number m < 0. This contrasts with the case of a Kerr black hole, which requires m > to 
realize > 0. This indicates the possibility of an unstable mode with m < 0. For m < 0, 
the condition P^ > is 

r+M 



-\mfl H \M 



qQ > fiM - qQ 



(34) 



for qQ > 0. However, the left-hand side of this inequality is negative, and thus Eq. 
cannot be compatible with the bound state condition (J2TJ). Thus for q > and m < in 
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Kerr-Newman spacetime, super radiance occurs, but the scalar field cannot be in a bound 
state, and the mode is stable. For Reissner-Nordstrom spacetime we have a = 0, and the 
relation > gives 

fi < — . (35) 
r+ 

This condition also cannot be compatible with the bound state condition (j21j) . and we 
conclude that there is no unstable mode of the scalar field in Reissner-Nordstrom spacetime. 
From these results, it is seen that the super radiance caused by the rotation of the black 
hole is essential to make the scalar field unstable. 

§3. Numerical approach 

3.1. Method 

To investigate the instability of the scalar field for a wide range of parameter values, we 
carried out numerical calculations. For this purpose, we first rewrite Eq. (J3J in terms of the 
tortoise coordinate 

.2 



X 



/r 1 
dr— = r H [r + ln(r — r+) — r_ ln(r — r_)] (36) 



and the new radial function u = rip. We then obtain 

d 2 u 

= Kff(r)n ' (37) 



Kff(r) 



A 



r 



2 



A 2M 2(a 2 + Q 2 ^ 



V 2 



1 

— [(r 2 + a 2 )uj — am — qQr]' 



A = /(/ + 1) — Imauj + [aio] 
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The effective potential V eS for a = 0.98M, Q = 0.01M, = 0.35, qQ = -0.08, I = m = 1 
is plotted in Fig-El Due to the mass of the scalar field, the effective potential has a well, 
and the wave can be trapped in this well. 

In the region near the horizon, the incoming solution of Eq. (J37)) is given by 



u ~ exp 



(38) 



j2Mu-qQ am + Q 2 cu\ 
V r + r l J 

and in the far region, the regular solution of Eq. (|37|) is given by 

,(M^_ ?Qa;)/v ^Z^ exp ^_ Xv / /i 2_ w 2j _ (39) 



We use the solutions (|38j) and (|3*^j) to impose the boundary conditions for the numerical 
integration of Eq. (|37j) . We prepare the inner numerical boundary x = X\ near the horizon 
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Fig. 2. The effective potential V e s(x) of the scalar field for a 
0.35, qQ = -0.08, I = m = 1. 



0.98M, Q = OMM, \iM 



r + and the outer numerical boundary x = x-i. By integrating Eq. (J37|) from x — X\ with the 
boundary conditions given by Eq. (JEEJ), we obtain a mode function vP"'. In the same way, 
we obtain with the boundary conditions imposed at the far region x = x-i- For a given 
complex value of uj, if the Wronskian 



W(u 



(i) u (2) : 



it 



(i) 



(2) 



(xi < x m < X2) is zero, the two solutions and u*- 2 -* 



(40) 



evaluated at the midpoint x = x m [Xi < x m < 22 J J s zero, trie two solutions u y "' and u y "> are 
linearly dependent, and u> is an eigenvalue of the equation under consideration. We search 
for the zero point of the complex function W(u) numerically in the complex uj plane. 

3.2. Result 

We performed numerical calculations to search for the mode of the scalar field satisfying 
I — m — 1 and uj ~ /i, for which the growth rate of the unstable mode is expected to 
have the largest value. We chose the parameter values of the black hole as a = 0.98M and 
Q = 0.01M and used the fourth-order Runge-Kutta integrator. The numerical boundaries 
were set at x\ = — 50M, x 2 = 1510M, x m = 24. 5M, and the grid spacing was Ax = 0.5M. 
We obtained the value of the growth rate 7 as a function of the scalar field mass fi and the 
scalar field charge q. We calculated the values of to for the eigenmode using 77 different sets 
of parameter values in g)-space. There points in (//, g)-space are indicated in FigO 

For the obtained value of to, we checked the validity of the assumption uj ~ \x used in 
Eq. (|37j) . In FigEl we plot \[i 2 /uj 2 — 1| as a function of \i for q = 0. For the obtained value 
of uj, this value is smaller than 5.5 x 10~ 2 , and our assumption uj ~ /i and (0) are correct 
for the numerically obtained modes. For q ^ 0, this value does not exceed 10~ 2 . We also 
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Fig. 3. The sets of parameter values used in the numerical calculations plotted corresponds to in 
(/i, g)-space. The solid line corresponds to = and the dotted line to M/j, — qQ = 0. The 
scalar field is expected to be unstable for sets of parameter values in the grey region. 
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Fig. 4. The value of |/U 2 /u; 2 — 1| as a function of \i for q = 0. This value is less than 5.5 x 10 -2 
for all values of /j,M, and we find that our numerical results are consistent with the assumption 

U> ~ jJL. 



calculated the eigen-modes using a different grid size. To check the numerical error on the 
obtained value of u, a numerical integration with x 2 = 3000M was also carried out. The 
relative error on the value of to evaluated from these calculations was found to be less than 
1(T 3 . 

The obtained growth rate is plotted in Fig.0 It is seen that the shape of the numerically 
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Fig. 5. The numerically obtained value of the growth rate 7. The spheres represent the obtained 
values of 7. The left panel is a birds-eye view of 7(/x, q) for a = 0.98M, Q = 0.01M, I = m = 1. 
The right panel is the same function as viewed from above. 



obtained function 7(/x, q) is almost the same as that of the analytically obtained one (see 
Fig.GJ. The left panel in Fig. |U] displays the /x dependence of 7 for different values of q. The 
plus symbols represent the maximum values of 7 for each /j, and the cross symbols represent 
7 for q = 0. The solid curve represents the values of 7 for q = predicted using the analytic 
approximation. In the region satisfying fiM < 0.25, the value of the numerically obtained 7 



0.8xlCr 7 













-0. 

-0.18 . 


-0.08 j.... -0.0. 


0.00 








■iQ- 




.22 £ 


















- 




t 


X 








-0.37 _^ 


- 




X 


























-0-55 

A 

++■-' 


'" 

> 















i.2xlo -7 
l.oxiir 7 

0.8X10 -7 
5 0.6xl0- 7 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 -0.7 











0.3 


o.:i2 .... 
) ,•**>■. 
+ + \ '»' 


0.35 

,(,.. 0.38 

1 V "$S 41 

V S ;•+ 








liM- 2 


0.28 .- 

*1 




\ t\ A * 
\ ' ' ,' \ ' 








; 


Yi :' 


+ 

y 


\ ?\ ' 










- 


20 / 

¥ \ 




; 


ft '+ 

I ...1.5, - 


fc - 
















H 


— *£ 


o.io A 

-■t' ++ ' / 








| \ - 


H 



-0.4 -0.3 -0.2 -0.1 0.0 0.1 



Fig. 6. The left panel displays the // dependence of the growth rate 7 for several values of q. The 
dashed curve represents the numerical solution for q = 0, and the solid curve represents the 
analytic approximation (|30|) for q = 0. The right panel plots the q dependence of the growth 
rate 7 for several values of [i. The dashed curves represent 7 for \i = 0.1 — 0.41. The solid 
curve represents the analytic approximation for fiM = 0.20. 
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approaches that of the analytic solution. The deviation from the analytic solution becomes 
significant for 0.3 < fiM. For parameter values in this region, the asymptotic matching 
method does not provide a good approximation. The right panel in Fig. El displays the q 
dependence of 7. The solid curve represents 7 obtained using the analytic approximation 
for fiM = 0.2, and it is seen that the approximation is good for —0.2 < qQ. 

The growth rate has a maximum value 7M ~ 1.13x 10~ 7 at fiM ~ 0.35, qQ ~ —0.08. The 
obtained minimum value of the growth rate is 7.6 x 10 -11 at fiM = 0.20, qQ = 0.1. Although 
the shape of the function 7(/i, q) agrees with the analytically obtained result displayed in 
Fig.^ its maximum value is three times larger. We confirm that the instability exists in 
the region of g)-space where both the super-radiant condition P^ > and the bound 
state condition Mfi > qQ are satisfied. This parameter region is also shown in Fig-El For 
all numerically obtained modes, the growth rates are positive, and they are contained in the 
region bounded by the two lines P^ = and Mfi — qQ = 0. As the parameter point (fi, q) 
approaches these lines, the growth rate decreases. We conclude that, the function 7(/i, q) has 
a maximum value in this region. We could not obtain a stable mode with negative 7; because 
the value of 7 for the stable mode is small compared to that for the unstable mode, it was 
not possible to obtain a definite value within the accuracy of our numerical calculation. 

In Fig.[7[ we show the behavior of the mode functions for the stable and unstable cases. 
The mode function of the stable mode m — — 1 increases monotonically while moving toward 
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Fig. 7. The behavior of the mode functions u(x). The solid curve is the real part and the dotted 
curve is the imaginary part of u{x). The left panel plots the mode function of the stable mode 
m = — 1. The right panel plots the mode function of the unstable mode m = 1. The parameter 
values are a = 0.98M, Q = 0.01M, fiM = 0.35, qQ = -0.08, 1 = 1. 



the horizon of the black hole. Contrastingly, the mode function of the unstable mode m = 1 
has a maximum at x ~ 27.5M. This location corresponds to the minimum of the effective 
potential (see Fig.|2J). For both values m — ±1, the effective potential has a well, which is 
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necessary in order for a bound state to exist. However, for the m = —1 mode, the super- 
radiant condition is not satisfied, and the wave falls into the black hole through the potential 
barrier. For the m — 1 mode, the super-radiant effect is significant, and the amplified wave 
can accumulate in the potential well. This leads to the instability of the mode. 

§4. Summary and discussion 

In this paper, we have studied the unstable modes of a massive scalar field in Kerr- 
Newman spacetime. We obtained the leading-oder value of the growth rate for \iM <C 
1, qQ <C 1 using the asymptotic matching method and for fiM < 1 using a numerical 
method. For a black hole with a = 0.98M, Q = 0.01M, we obtained the maximum value 
of the growth rate of the unstable mode to be 7M ~ 1.13 x 10 -7 for fiM ~ 0.35, qQ ~ 
—0.08. The location of the maximum value in g)-space agrees with the result of the 
analytic method, but its numerical value is three times larger than that of the analytic 
result. For 0.3 < /iM, the numerically obtained value of 7 deviates significantly from that 
of the analytic result . This indicates that the asymptotic matching method does not give a 
good approximation in this parameter region. If the mass of the scalar field is fi ^> l[MeV], 
the mass of the black hole must be M < 10 15 [Kg] for that the instability to be maximal. 
Thus, there is a possibility that the mechanism of the black hole bomb can influence the 
evolution of primordial black holes. 

To investigate the dynamics of the unstable mode, it is interesting to consider the tem- 
poral evolution of a scalar wave with an unstable mode in the black hole geometry. The 
propagation of such a wave in black hole spacetime consists of three stages. During the first 
stage, there is a burst wave that depends on the initial conditions of the wave. This stage is 
followed by quasi-normal ringing and the tail mode. If the scalar field has an unstable mode, 
we expect the effect of the instability to appear in the late time tail behavior. However, 
because the growth rate of the instability is very small, it may be difficult to detect the 
instability numerically. This is the next problem we intend to study. 
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